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abstract 

m ' 

Context. The star HD 49385 is the first G-type solar-like pulsator observed in the seismology field of the space telescope 
CoRoT. The satellite collected 137 days of high-precision photometric data on this star, confirming that it presents 
solar-like oscillations. HD 49385 was also observed in spectroscopy with the NARVAL spectrograph in January 2009. 
^ , Aims. Our goal is to characterize HD 49385 using both spectroscopic and seismic data. 

5^ ■ Methods. The fundamental stellar parameters of HD 49385 are derived with the semi-automatic software VWA, and 

the projected rotational velocity is estimated by fitting synthetic profiles to isolated lines in the observed spectrum. 
A maximum likelihood estimation is used to determine the parameters of the observed p modes. We perform a global 
fit, in which modes are fitted simultaneously over nine radial orders, with degrees ranging from t — to I — 3 (36 
individual modes). 

Results. Precise estimates of the atmospheric parameters (T e ff, [M/H], log<?) and of the vsini of HD 49385 are obtained. 
The seismic analysis of the star leads to a clear identification of the modes for degrees I = 0, 1, 2. Around the maximum 
of the signal (y ~ 1013 ^Hz), some peaks are found significant and compatible with the expected characteristics of I = 3 
■ modes. Our fit yields robust estimates of the frequencies, linewidths and amplitudes of the modes. We find amplitudes 

of ~ 5.6 ± 0.8 ppm for radial modes at the maximum of the signal. The lifetimes of the modes range from one day 
(at high frequency) to a bit more than two days (at low frequency). Significant peaks are found outside the identified 
, ridges and are fitted. They are attributed to mixed modes. 

Key words. Methods: data analysis - Methods: statistical - Methods: observational - Stars: oscillations - Stars: indi- 
vidual: HD 49385 

1. Introduction Sun, allowing us to estimate t he sound speed and den - 

' T .„ . . , , , sity profiles (iBasu et al. 20031 iTurck-Chieze et al. 2 001). 

H ; In the Sun, escalations are excited by the turbu- the position of the base of the convective zone 

lent motions m the outer part of the external con- ( Christcnscn-Dal sgaard et al. 1991. and the ro tation pro- 

vective envelope and are further propagated into the file (: T hom P s oD^t al. 2003L iMathur et al. 20081) . However, 

interior of the star. The study of these oscillations the yery bw amplitude of these oscil i ations ( a f ew ppm m 

has yielded constraints on the inner structure of the photometr y) makes it very challenging to detect and ana- 

~ Z ~ ~~ , , „ „ , , lyze them in other stars than the Sun. 
bend offprint requests to: S. Deheuvels 

e-mail: sebastien.deheuvels@obspm.fr a i • u j , j- c ,i • , r 

. r , , . Achieving a better understanding oi the interiors oi 

Based on data obtained from the CoRol (Convection, , , ,. c . , ,. r , , 

„ , , . , . _ . . \ . . ) . , ' solar- like pulsations is one of the mam obiectives ol the 

Rotation and planetary transits) space mission, developed by . . „ _ _ ,. _ ; . . 

the French Space agency CNES in collaboration with the Science s P ace missi0n CoRoT (Convection, Rotation and plan- 
Programs of ESA, Austria, Belgium, Brazil, Germany and etai T Transits). CoRoT is a space telescope performing 
Spain. high-precision p hotometry over qu asi-uninterrupted long 
** Based on data obtained using the Telescope Bernard Lyot observing runs (|Baglin et al. 2006|). Solar-like oscillations 
at Observatoire du Pic du Midi, CNRS and Universite Paul have already been studied in several other stars with 
Sabatier, France. CoRoT data. The analyses of these stars have encoun- 
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tered difficulties identifying the degrees of the modes, ei- 
ther because of a to o low signal- to-noise ratio (HP 175726: 
iMosser et al. 2009U HD 181906: iGarcia et al. 20091 ). or be- 
cause of a too short lifetime of the modes, inducing large 
mode line widths (HD 49933: lAppourchaux et al. 20081 
HD 181420: iBaTban et al. 20091 ). 

The star HD 49385 is the first G-type solar-like pulsator 
observed in the seismology held of CoRoT. It is cooler than 
the solar-like pulsators previously analysed with CoRoT 
data, and probably more evolved (at the end of the Main 
Sequence or shortly after it). The choice of HD 49385 as 
a CoRoT target has motivated us to lead spectroscopic 
observations, performed with the NARVAL spectrograph 
at the Pic du Midi Observatory. The fundamental param- 
eters of HD 49385 are derived from these observations, as 
described in Sect. The photometric observations with 
CoRoT are presented in Sect. [31 Section @] presents the 
study of the low-frequency part of the power spectrum, in 
search of a signature of the stellar rotation. The extraction 
of p-mode parameters is described in Sect. [SI and Sect. [5] is 
dedicated to conclusions. 



2. Stellar parameters of HD 49385 

The solar-like pulsator HD 49385 is a GO-type star with 
an apparent magnitude o f my = 7.39 (uvby catalog, 
lHauck fe Mermilliod. 1 998) . To determine the fundamental 
parameters of HD 49385 we analysed a high-quality spec- 
trum from the NARVAL spectrograph mounted on the 2- 
m Bernard Lyot Telescope at the Pic du Midi Observatory. 
We stacked two spectra collected on 2009 January 10, total- 
izing an exposure time of 6 600 s. The spectrum was normal- 
ized by identifying continuum windows in a synthetic spec- 
trum and fitting a low-order spline through these points. 
This was done order-by-order and we made sure the line 
depths agreed for the overlapping part of adjacent echelle 
orders. We measured a mean signal-to-noise (S/N) ratio of 
600 in the continuum for several line-free regions in the 
range 5 000-7 000 A for a data sampling of 2.7 data points 
per resolution element (i? = 65 000). 



2.1. Temperature, log g and metalicity 

We used the semi-automatic VWA software ([Bruntt 2 009) 
to fit synthetic profiles for more than 600 lines in the 
range from 4135 to 8 545 A. We used a differential ap- 
proach, meaning all abundances are measured relative 
to e xactly the same li nes in the solar spectral atlas 
from lKurucz et al. (19841 ). This differential ap proach is de- 
scribed in more detail by iBruntt et al. (20081 ) . As part of 
the analysis, we determined the atmospheric model param- 
eters T e ff, log<7, and the microturbulence. This was done 
by requiring that abundances determined from Fe I lines do 
not correlate with the equivalent width (EW) or the ex- 
citation potential (EP). Furthermore we required that the 
same mean abundances are measured from the Fe I and Fe II 
lines ("ionization balance"). We estimated the uncertainty 
on the model parameters by adjusting them until the cor- 
relations of Fe I with EW or EP became significant or the 
ionization balance deviated (see IBruntt et al. 20081 for de- 
tails). We stress that the uncertainties are strictly internal 
errors because the underlying assumption is that the model 
atmosphere represents the star. Any systematic error in the 



Table 1. Spectroscopic parameters of HD 49385. 



T cB 6095 ± 65 K 

log g 4.00 ±0.06 

[M/H] +0.09 ± 0.05 dex 

vsmi 2.9±i;| 



kms" 



temperature profile or departures from local thermal equi- 
librium (LTE) will likely affect the results. We therefore 
quadratically added 50 K and 0.05 dex to the uncertainty 
on T e g and log g. The spectroscopic parameters are listed 
in Table [U 

As an additional check of the surface gravity we fit- 
ted the wings of the pressure sensitive Ca lines at 6122.2 
and 6162.2 A. This was done by first adjusting the van der 
Waals constants to fit the solar spectrum for the canonical 
value \ogg — 4.438. The fitted values for the two lines for 
HD 49385 are log 3 = 3.96 ± 0.12 and 4.05 ± 0.06, which 
agrees very well with the Fei/Fen ionization balance. The 
low value we obtain for the log g of HD 49385 indicates that 
it certainly is an evolved object. 

The parameters agree well with the c alibration of 
the Stromgren indices From TEMPLOGG ([Rogers 1995L 
iKupka fe Bruntt 20011) we get T eff = 6345 ± 150 K, log 5 = 
4.22 ± 0.15, and [Fe/H] = +0.21 ± 0.10. TEMPLOGG 
determines a signihcant interstellar reddening E(b — y) = 
0.020 ± 0.007. If we assume z ero interstellar reddening, we 
get from the V — K, using the[R amircz fe Melendez (2 005) 
calibration, T cS = 5958 ± 64 K and 6035 ± 92 K from 
b — y. The uncertainties include the calibration uncertainty 
and the photometric uncertainty on the indices. Finally, 
iKovtvukh et al. (2004h used line-depth ratios to determine 
effective temperatures of stars and found T e g = 6052+6. 7K 
(internal error). They also determined log g and [Fe/H], 
but did not give uncertainties on these two parameters: 
log g — 4.00 and [Fc/H] = +0.10. In summary, our result 
from VWA agrees well with the photometric calibrations 
and the line-depth ratio method for T e g. 

We found a high abundance of lithium in HD 49385 
of +1.08 ± 0.10 (relative to the Sun). This value ap- 
pears to b e typical for a dwarf star with solar metal- 



icity Je^_ Boesgaa^^r^r^OO^) . We used the line list 



from lGhezzi et al. (20091) . but did not include the relatively 
weak molecular CN bands; hence the Li abundance may 
be slightly overestimated. The Li feature around 6707.8 is 
known to be strongly affected by NLTE effects, but when 
using 3D instead of ID atmospher e models these two effects 
appear to cancel to first order (|Asplund 20051) . In other 
words, because we are using ID LTE models, the imposed 
error should be relatively small. Furthermore, because we 
are calculating abundances relative to the same Li line fea- 
ture as in the Sun, our estimate of [Li] should be robust 
as the atmospheric parameters of the two stars are quite 
similar. 

The abundance pattern we determine is shown in Fig. [T] 
and listed in Table O The overall metalicity is determined 
as the mean of the metals with at least 10 lines used in 
the analysis (Si, Ca, Ti, V, Cr, Fe, Co, Ni) giving [M/H] = 
0.09 ± 0.05. This range is marked by the horizontal bar in 

Fig.m 
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Fig. 1. Relative abundances of 22 elements measured in 
HD 49385. Open and solid circles are mean values for neu- 
tral and singly ionized lines, respectively. The yellow hori- 
zontal bar marks the mean metalicity with 1-a uncertainty 
range. 



2.2. Projected rotational velocity 

Because we have a relatively high-resolution spectrum with 
very high S/N we were able to estimate the projected 
rotational velocity (vsini). We fitted synthetic profiles 
to ten isolated lines in the range 4 600-6 800 A. We cal- 
culated the synthetic pro files with the SYNTH software 
IValenti fe Piskunov (199(f ). thus taking into account weak 
blends. We considered three broadening mechanisms which 
were convolved with the synthetic profiles: (1) the instru- 
mental resolution, (2) macroturbulence (ii macro ) and (3) 
v sin i. The instrumental power of resolution of NARVAL is 
R = 65 000, but to check this value we fitted six narrow tel- 
luric lines around 6290 A. We found that a Gaussian func- 
tion with a FWHM of 6290 A/65 000 = 0.097 A fitted these 
lines very well and we used the same value of R for the stel- 
lar lines. For -usini and i> maC ro we calculated a 17 x 17 grid 
of convolved profiles for 0-8 kms -1 in steps of 0.5 kms -1 . 
To identify the best combination of v sin i and v ma cro we 
compared the synthetic profile and the observed spectrum 
by computing the x 2 value. 

In Fig. [2] we show an example of the \ 2 contour for 
the Fei line at 6156.1 A. The plots for all ten lines are very 
similar and show the same strong correlation between v sin i 
and Wmacro- We marked three points on the contour in Fig. [5] 
with filled circles, and the synthetic profiles corresponding 
to these grid points are shown in Fig. [3] The observed spec- 
trum is shown with open circles. From our analysis we can 
put a firm upper limit on vsini of 5 kms -1 , corresponding 
to Umacro = kms -1 . The values that provide the best fit 
based on the analysis of the \ 2 contours for all ten lines are 
usini = 2.9^'g kms -1 and w ma cro = 3.4^q'y kms -1 . 

2.3. Luminosity 




2 4 6 

v sin i [km/s] 



Fig. 2. Example of the \ 2 contour surface used to deter- 
mine v sinz and u m acro for the Fe I line A6151.6 A. The three 
circles mark the minimum of the surface and the 1-a un- 
certainty range. The corresponding three synthetic spectra 
are compared to the observed line in Fig. [31 




-0.02 r i 
-0.03 1 . , . . . , . . . , . . . , . . . I 

6151.2 6151.6 6152.0 

Wavelength [A] 



Fig. 3. The upper panel shows the fit of three synthetic 
profiles to the observed line (marked by open circles). Each 
profile was convolved with different combinations of v sin i 
and i> maC ro as given in the top panel (cf. Fig. [5]) . The lower 
panel shows the difference between the observed and syn- 
thetic profiles, indicating a slightly better fit for the profile 
represented by a dashed line (the \ 2 value for each fit is 
given in the lower panel). 



Usin g the Hipparcos par allax of the star, ir = 13.91 ± 0.76 
mas ()van Leeuwen 20071 ) . we could derive the absolute vi- 
sual magnitude of the object My = 3.11 ±0.12. The 
bolometric correction was determined by interpolating in 
the grid provided by iBessell et al. (19981 ) for appropriate 
values of temperature and \ogg for HD 49385. We ob- 
tained BC\ — —0.029 ± 0.006, which yields a luminosity 
log(L/Z ) = 0.67 ± 0.05. The values of T cff and L provided 
an estimate of the radius R/Rq = 1.94 ± 0.15. 



3. CoRoT observations 

The star HD 49385 was one of the targets of the second 
Long Run in the asteroseismology field of the space mission 
CoRoT. 136.9 days of photometric data have been col- 
lected from October 2007 to March 2008, with a duty cycle 
of 88.2%. After bein g corrected from know n instrumental 
effects as described in lSamadi et al. C2007al ). the lightcurve 
consists of a series of equally-spaced measurements (32 s) in 
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Table 2. Abundances of 22 elements in HD 49385 relative 
to the Sun. The third column lists the number of lines used 
to determine the mean abundance listed in the second col- 
umn. The quoted uncertainties are intrinsic RMS errors on 
the mean value. 



El. 




n 


El. 


AA 


n 


C i 


-0.03 ±0.05 


9 


V i 


+0.09 + 0.03 


12 


N i 


+0.18 


1 


Cr I 


+0.09 ± 0.03 


18 


O i 


+0.18 ±0.17 


4 


Cr II 


+0.08 ± 0.03 


7 


Na I 


+0.08 ± 0.04 


4 


Fe I 


±0.11 ± 0.03 


336 


Mg i 


±0.08 ±0.14 


4 


Fe II 


±0.11 ±0.03 


27 


Mg n 


+0.07 + 0.03 


2 


Co I 


±0.03 ± 0.03 


12 


Al i 


+0.10 + 0.03 


4 


Ni i 


±0.06 ± 0.03 


81 


Si i 


+0.08 ± 0.03 


38 


Zn I 


-0.04 


2 


Si ii 


±0.06 ± 0.04 


2 


Sr I 


±0.09 


1 


S i 


-0.02 


2 


Sr II 


±0.10 


1 


Ca i 


+0.14 + 0.03 


13 


Y II 


±0.10 ±0.04 


7 


Ca ii 


+0.16 + 0.03 


2 


Ba ii 


±0.27 


1 


Sc ii 


+0.15 + 0.04 


6 


Ce ii 


±0.06 


1 


Ti i 


+0.08 + 0.03 


30 


Li i 


±1.08 


1 


Ti II 


+0.17 + 0.03 


9 









the heliocentric frame (also known as N2 data). Small gaps 
remain in the lightcurve (a few minutes long) mainly due to 
the passage of the CoRoT satellite in the South Atlantic 
Anomaly (SAA). A larger gap (3.5 days) was caused by 
an unexpected reset of the DPU (Digital Processing Unit), 
which was probably the consequence of the impact of an en- 
ergetic particle, because this event occured while the satel- 
lite was crossing the SAA. These gaps were linearly interpo- 
lated, based on the measurements collected 1 000 s on either 
sides of the gap. A low-frequency trend due to the aging 
of the CCD was also corrected b y removing a linear g ain 
variation of —5.32 10 -5 day -1 (see lAuvergne et al. 20091 for 
more details). The resulting lightcurve is shown in Fig. 01 
Figure [5] shows the power spectrum of the lightcurve, 
computed with the Fast Fourier Transform algorithm 
(FFT). It was normalized so that the integrated power 
spectrum from to twice the Nyquist frequency corre- 
sponds to the variance of the time series. At low frequency 
the background rises due to stellar granulation and possi- 
bly to stellar activity (see discussion in Sect. 01 . The sig- 
nature of perturbations due to the orbit of the satellite 
(v 01 h — 161.7/iHz) and the day (vday — 11.6 /xHz) remains 
in t he power spectrum. Th ese perturbations were expected 
(see lAuvergne et al.. 20091 ) and are of relatively low ampli- 
tudes (a few tens of ppm 2 /^Hz for the strongest ones at 
low frequency). They appear as combinations of the har- 
monics of both frequencies. For instance, a peak appears at 
150.1 /xHz (corresponding to ^ or b — t'day)- Finally, a broad 
excess of power, centered around 1 mHz, is clearly seen in 
the power spectrum and corresponds to the acoustic modes 
of HD 49385, which are studied in Sect. [SJ 



4. Search for a signature of stellar rotation in the 
Fourier domain 

The measurement of the rotation period from the lightcurve 
analysis was found to be us eful in previous analy- 
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Fig. 4. Lightcurve of 137 days of CoRoT observations on 
HD 49385, detrended using instrumental information. 
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ses of p-mode par a meters (see Appourchaux et al. 20081 



iBarban et al. 2009L iMosser et al. 2009bl) . It indeed pro 



Fig. 5. Power spectrum of 137 days of CoRoT data for 
HD 49385. The peaks indicated by arrows are from left 
to right: the 14-^Hz-peak discussed in Sect. HI a peak at 
2^day = 23.1 /^Hz, and a peak at i/ or b — Vday = 150.1 /xHz. 
The grey area shows the excess of power due to the p-mode 
oscillations. 



vides a direct measurement of the surface rotation period 
and hence a first estimate of the rotational splitting value. 

Contrary to other solar-like pulsating CoRoT targets, 
no immediate signature of the activity can be seen in the 
lightcurve (see Fig. 2]). Only a long period trend stands out. 
The same trend was also found in several stars observed on 
the same CCD during the same specific observation run. 
Its instrumental or environmental origin is thus established 
and is currently under study. 

We first identified a peak around 0.4 fiRz (period of 
about 29 days) , the significance of which is hard to assess 
because it presents a relatively low amplitude in a frequency 
domain where the intrumental noise can dominate the 
power spectrum. We performed a tim e-frequency analysis 
with Morlet wavelets (as described in lMathur et al. 2 010) 
which showed that this peak presents a certain stability in 
time. But because the periods we were testing correspond 
to a non-negligible fra ction of the total observat ion period, 
the cone of influence (jTorrence fe Compo 19981) prevented 
us from drawing any conclusion. 

A prominent peak can also be detected at a frequency 
of 14 yitHz (period of about 0.8 day) with a signal-to-noise 
ratio of more than 10. This peak is however unlikely to 
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400 600 800 1000 1200 1400 1600 
Frequency (/iHz) 

Fig. 6. Smoothed power spectrum of HD 49385 (with a 10 
/iHz boxcar) in the frequency range of the oscillations. The 
green curve represents the fitted background (granulation 
and white noise components) . The blue dashed curve corre- 
sponds to a Gaussian fit of the raw power spectrum giving 
^max = 1013 /iHz (dotted line). 



be directly related to rotation. Indeed, let us assume ad 
absurdum that the rotation frequency is 14 /iHz. Our esti- 
mate of the radius implies a surface velocity of more than 
120 kms -1 . To obtain the measured -usini, we would need 
a value of the inclination angle lower than 2°, which means 
that the star is seen pole-on. In this case it is extremely 
difficult to derive the signature of the activity from the 
light curve (M osser et al. 2009ah and no peak related to 
the rotation period can be detected in the power spectrum, 
which contradicts our assumption. In the power spectrum of 
HD 49933, which is a target of the seismo- field for the same 
run, we also found a small peak around 14 /iHz, suggesting 
this peak might be instrumental. However, the frequencies 
do not match exactly, and the peak observed for HD 49933 
does not appear to be significant. The origin of this peak 
remains unclear. But it cannot be related to rotation. 

In summary, we found no clear signature of the stel- 
lar rotation of HD 49385. The low- frequency part of the 
power spectrum can therefore not be used to obtain a first 
estimate of the rotational splitting. We may however infer 
from the absence of a clear activity signal either that the 
star has a low activity, which cannot be detected even with 
the CoRoT sensitivity, or that the inclination of the star 
is close to pole-on. 

5. Analysis of p mode oscillations 

Figure [5] shows the power spectrum smoothed with a 10- 
/iHz boxcar. The comb- like structure of this excess of power 
is typical of p-mode oscillations. This section presents the 
analysis of these modes. 

5.1. Frequency of maximum power and large spacing 

To estimate the frequency of the maximum of signal, we 
fitted a Gaussian profile on the raw power spectrum in the 
frequency range of the oscillations and above a fitted back- 
ground (fit described in Sect. I5.3|) . The result is shown in 
Fig. [B] We obtained as the frequency of maximum power 

'-'max — 1 013 ± 3 /iHz. 



n 




Frequency modulo Al/ (/iHz) 



Fig. 7. Echelle diagram of HD 49385 computed with a large 
spacing of = 56 /iHz. The power spectrum is smoothed 
with a l-/iHz boxcar and normalized by the background (es- 
timated in Sect. I5.3.T]) . Three clear ridges appear, labeled 
here as ridges A, B, and C. The white dashed lines corre- 
spond to the frequencies obtained from a fit to the data 
presented in Sect. [5751 and summarized in Tables [3] and |H 



In Fig. [5] a comb-like structure is clearly visible between 
600 and 1400 /iHz. The autocorrelation of the power spec- 
trum provides a first estimate of the mean value of the large 
spacing Av ~ 56 /iHz. 

5.2. Mode identification 

We used the value of the large spacing given in the previous 
section to build the echelle diagram shown in Fig. [7J It is 
obtained by piling onto one another sections of 56 /iHz of 
the power spectrum. Three clear ridges appear in the echelle 
diagram: two neighboring ones, referred to as ridges A and 
B in Fig. [JJ and a third one (ridge C). 

The pattern we observe in Fig.[7Jis the one we expect to 
observe for low-degree high-radial-order modes in solar-like 
pulsators: two main ridges separated by about half a large 
spacing corresponding to £ = and I = 1 modes, and 
alongside the I — ridge a fainter one corresponding to 
I — 2 modes. It then seems reasonable to identify the 
ridges as follows: neighboring ridges A and B correspond 
to i = 2 and t = mode s, respectively, and ridge C to 
I = 1 modes. We note that | Mosser fc Appourchaux (2 009) 
and iBedding fc Kieldsen (201(1 " obtained the same iden- 
tification for the degrees of the ridges using different 
approaches. The following analysis is based on this identi- 
fication. However, since the previous analyses of CoRoT 
solar-like tar gets encountered ambiguity indentifyin g the 
modes (see iBarban et al. 20091 iGarcia et al. 2009 i. we 
decided to also consider the possibility of the alternate 
identification (for which ridge C corresponds to I = 
modes), which is discussed in Sect. 15.3.41 

Figure [Jj also shows fainter peaks outside the identi- 
fied ridges. We here study their significance by testing 
the Ho hypothesis in the frequency range of the oscil- 
lations. Because we searched for short-lived modes, we 
binned the power sp ectrum over n bins, as prescribed in 
App ourchaux f2004tl . The size of the boxcar (n = 15 bins 
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~ 1 /iHz) for the binning was chosen to approximately 
match the expected linewidth for the p modes in solar- 
like pulsators. In this case, the noise is distributed as a \ 2 
with 2n degrees of freedom. We then applied the Ho hy- 
pothesis, i.e. we supposed that the observations are due to 
pure noise. For each peak we computed the probability p 
that pure noise is able to reach the height of the peak over 
the whole studied window. Low values of p indicate a poor 
compatibility between the observations and noise, and the 
significance of the peaks is defined as 1—p. We show in Fig. 
[5] all the peaks with a significance higher than 95%. 

We notice four peaks which appear with a significance 
above 98%, and lie outside the identified ridges. They are 
labeled as ttl, 7r 2 , 7r 3 and 7r 4 on Fig. [5] The peak 7T4, clearly 
outside the pattern expected from asymptotic theory, could 
correspond to a mixed mode, given that HD 49385 is an 
evolved object. The peak 7Ti might be an I = mode, as- 
suming that the curvature of ridge B changes at low fre- 
quency. This hypothesis seems to be corroborated by the 
fact that such a change exists in ridge C. However, a more 
precise inspection shows that including ix\ in ridge B gener- 
ates an abrupt step in the £ = large separation, whereas 
the £ — 1 large separation varies quite smoothly. The peak 
7Ti could also be an I — 1 mixed mode in avoided crossing. 
Indeed, at low frequency, the £ — and 1=1 ridges have 
different curvatures. This can be seen in Fig. [7] and will 
be confirmed by the large separation profiles derived from 
th e fitted mode frequencies (s ee Sect. I5.3."2l) . It was shown 
in iDeheuvels fc Michel (20091) that a low-degree i avoided 
crossing creates a characteristic distortion in the ridge of 
degree £. If we assume that tt\ is an £ = 1 mixed mode 
in avoided crossing, the expected curvature of the £ = 1 
ridge is very similar to that of the observed ridge. We can- 
not establish this identification more firmly at this stage 
though. Because of the uncertainty regarding the identifi- 
cation of the 7Ti mode, we preferred to consider and fit it 
individually (see Sect. I5.3[) . The peaks 7r 2 and n 3 are lo- 
cated at about 7 /iHz on the left of ridge C. We studied the 
possibility for these peaks to correspond to mixed modes. 
In this case, they would have to be of different degree £. 
Indeed, g modes of same degree and increasing radial or- 
der have frequencies increasingly close to each other. Two 
mixed modes of same degree £ spaced by Av would imply 
many more mixed modes in the ridge of degree £ at lower 
frequency. It is clear from ridges I — 1 and I — 2 that this is 
not the case. And if 7r 2 and 7T3 were mixed modes of different 
degree, we would expect them to have different amplitudes, 
which does not seem to be the case. Their position in the 
echelle diagram on the left of ridge C suggests that they 
could also be £ = 3 modes. Several reasons led us to favor 
this identification. First, these two modes appear vertically 
aligned in the echelle diagram and around the maximum 
of the oscillations, i.e. at the place where it is most likely 
to observe £ — 3 modes if they are present. Then, 7r 2 and 
7T3 might be part of a fainter ridge (referred to as ridge D), 
because other peaks, which have a lower significance but 
seem to follow the alignment, appear in Fig. [7J Finally, the 
distance between ridges C and D is also consistent with 
this identification. Indeed, the asymptotic theory gives an 
expression for the average distance between the ridges 

(6is n . e ) = v n< i - i/ n _ ljt+2 = (4£ + 6)D , (1) 

where the expression of Dq is found in 
IChristens cn-Dalsgaard ~(2003l ). By definition, the quantity 
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Fig. 8. Echelle diagram of the peaks with a significance 
greater than 95% in the power spectrum smoothed over 
n = 15 bins (~ 1 /iHz). Four peaks appear, which do not 
belong to the identified ridges A, B, or C, and are labeled 

as 7Ti, 2 ,3,4. 

(Sv nt e=o) corresponds to the distance between the I = 
and £ — 2 ridges, and {8v n ,i=i) corresponds to the distance 
between the £ = 1 and £ = 3 ridges. Figure [7] shows that we 
have (5i>n,e=o) — 5 /iHz. We therefore expect I — 3 modes 
to be approximately 8 /iHz on the left of £ = 1 modes, 
which approximately matches the position of ridge D. 

The significance is an efficient tool to select peaks, but 
it is not a guarantee for the d etection of a signal, as stressed 
in lAppourchaux et al. (20091) . To confirm the existence of 
a signal, one needs to compute the posterior probability 
P(Hq\X = x). The details of the calculation along with 
the assumptions we make are given in Appendix [Bj For 
iv i and 7T4, assumed to be mixed modes, we found poste- 
rior probabilities of 0.001% and 8.7%, respectively. These 
peaks are therefore probably due to signal. The larger pos- 
terior probability found for 7T4 is caused by the fact that the 
expected profile of mixed modes was built with very con- 
servative priors . For 7r 2 and tt^, assumed to correspond to 
£ — 3 modes, we obtained a very low posterior probability 
(1.4% and 1.1%, respectively), indicating that these peaks 
are highly likely due to signal, and a good compatibility 
between their observed profiles and the expected profile of 
I — 3 modes for HD 49385. This confirms our identification 
for 7r 2 and 773. It is the first time that £ = 3 modes can 
be detected in a solar-like pulsator (other than the Sun). 
We explain in Sect. 15.3^21 why HD 49385 was a particularly 
favorable target to detect £ = 3 modes. 

We remark that ridge D could also be seen as the m 
■ 1 component of the rotationally splitted £ = 1 mode. This 
hypothesis is studied in 15.3.41 

5.3. Extraction of p mode parameters 
5.3.1. Method 

Solar-like oscillations are the result of stochastic excita- 
tions of p modes in the convective envelope. They re- 
spond to the equation of a forced , damped and randoml y 
excitated harmonic oscillator (see lDuvall fc Harvev 1 986). 
Therefore their signature in the power spectrum can 
be modeled as Lorentzian profiles with a linewidth pro- 
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portional to the inverse of the mode lifetimes, added 
to the background and perturbated by a multiplica- 
tive noise following a x 2 distribution with two degrees 
of freedom. We used a maximum likelihood estima- 
tion (MLE) to determine the m ost probable parame - 
ters of the observed p modes (see I Anderson et al. 1 990). 
This method was already applied to determine p- 
mode parameters for the Sun ([Toutain fe Froehlich 1 992. 
Appourc haux et al. 19981) as well as for CoRoT solar- 
like targets (lAppourchaux et al. 20081 Barban et al. 2009. 



value of the limb darkening for this star into account, we 



iGarcia et al. 2009(1 . As described in Anderson et al. f 1 9901 , 
the likelihood L is defined as the product of the probability 
density of each frequency bin vi of the power spectrum 



L(v l ,\) = Yl 



exp 



-S(vi 



M{v u X) 



(2) 



where S is the observed power spectrum, M the model 
described above, and A the parameters of the Lorentzian 
profiles used in the model. The set of parameters A that 
maximizes the likelihood corresponds to the most probable 
parameters. Maximizing L is the same as minimizing the 
negative logarithmic likelihood C, defined as 



C = -ln(L) = (hiM(i/i,X 



M{yi,X) 



(3) 



As mentioned above, perturbations remain in the power 
spectrum at frequencies corresponding to combinations of 
the harmonics of the orbital frequency and those of the day. 
Some of these perturbations can be seen in the frequency 
domain of the oscillations. We therefore gave a null weight 
to the corresponding bins so that our estimates of the p- 
mode parameters are not influenced by those peaks. 

We estimated the background with a power law to 
describe the granulation component, as prescribed by 
lHarvev (1 985). Prior to the fit of the p-mode parameters, 
we fitted a profile of the type 



B(ui) = 



1 + (/^) 7 



(4) 



on the power spectrum, excluding the frequency domain of 
the p-mode oscillations to avoid any bias, with a, /?, 7 and 
S as free parameters. We then used this background to fit 
the mode parameters. 

The analysis of HD 49933, which was the first solar-like 
CoRoT target, showed that a global fit of all the modes 
gave more robust estima tes of the mode parameters (see 
Appourc haux et al. 20081 ). We therefore fitted modes simul- 
taneously on a frequency range corresponding to nine over- 
tones in the power s pectrum. The assumptions m ade for the 
fit are the same as in lAppourchaux et al. (200"8j) . We briefly 
recall them here. The splitting was assumed to be con- 
stant over the frequency range of the oscillations, and the 
m 7^ components of non-radial modes were supposed to 
be symmetrical with respect to the m = component. We 
sliced the power spectrum in A^-wide intervals, in which 
the linewidth is assumed to be the same for all modes, and 
where only the height of the £ — mode (-H^=o,n) is fitted. 
The heights of non-radial modes are determined from those 
of the £ = mode. Assuming an equal repartition of the 
energy in the modes of different degrees £, and taking the 



obtain 




H-l=i,nl H-t=v,n 


= 1.5 


Hg=2,n/ He=o,n 


= 0.5 




= 0.05, 


where Hi—\^ n 





>H n j_ m . Inside non-radial multi- 
plets, th e heights H n .i.m are given by the intensity visibili- 
ties (see lGizon fc Solanki 20031 iBallot et al. 200(4 



5.3.2. Mode parameters estimates 

Five teams worked independently on the extraction of the 
p-mode parameters and applied the method we described 
in Sect. l5.3Tl to fit the four ridges associated to I = 0, 1, 2, 3 
modes. As was done in previous analyses of this type, we 
present in the following section one set of reference results, 
and we use results from the other teams to cross-check them 
and comment the robustness of the solution. The reference 
results of the fit are given in Tables and S) 



Mode frequencies 

There is a very good agreement between the different 
teams on the fitted mode eigenfrequencies. Only two 
frequencies in the t = 3 ridge differ by more than 1 a. 
Indeed, this ridge is too faint to offer reliable estimates 
of the mode frequencies away from the maximum of the 
signal. An echelle diagram of the fitted modes is given 
in Fig. |TD] At the maximum of the signal, the obtained 
precision on the mode frequencies ranges from about 0.2 
to about 0.6 /iHz for ^ I ^ 3 modes and increases at the 
edges of the frequency domain of the oscillations. 

At v = 855.3 /iHz, the 1 = 1 mode is found merged 
together with the I = mode. This could be because the 
signal-to-noise ratio decreases when reaching the edges of 
the fitted domain. Indeed, with a lower signal-to-noise ra- 
tio, the analysis of previous CoRoT targets yielded sev- 
eral £ = 2 modes merged with th e neighboring £ = 
mode (see lAppourchaux et al. 20081 iBarban et al. 2 009). 
However, the £ — 2 mode might also have been shifted 
because of an £ = 2 mixed mode in the neighborhood. 

Figure[5]shows the large separation profiles derived from 
the fitted frequencies of the £ = and £ = 1 modes. It 
appears, as mentioned in Sect. 15.21 that the £ = and 
£ = 1 ridges have a different curvature at low frequency, 
which might be caused by a low-frequency avoided crossing. 

Mode linewidths 



Figure [Til shows the linewidths of the fitted modes. As 
was observed in the Sun, the linewidths of the modes in- 
crease with frequency and are almost constant around the 
maximum of the signal. They are larger than those of the 
Sun by a factor 2. The lifetimes r of the modes can be de- 
duced from the linewidth by the relation r = l/(7rr). For 
HD 49385 they range from about one day to two days. By 
comparison the mode lifetimes are shorter in other CoRoT 
solar-like pulsators, e.g. for HD 181420 (about six hours for 
the shortest ones, and 12 hours around the maximum of 
the signal). 
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Table 3. Reference results for the parameters of the radial p modes. The heights given here are those of the £ = modes 
(the heights of non-radial modes are computed using height ratios, see Sect. I5.3TT]) . The agreement (agr.) between the 
teams is given for each fitted parameter. Filled circles indicate that at least four of the five teams agreed on the value 
within 1 a. X a indicates that more than one team disagrees with the value, and that the five teams agree within X a 
(X > 1). 
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Fig. 9. Large separation profiles derived from the fitted fre- 
quencies of the I = (squares and blue line) and I = 1 
(triangles and red line) modes. We overplotted the 1-a er- 
ror bars. The cross corresponds to the large separation we 
would obtain, if wi were considered as an £ = mode. 
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Fig. 10. Echelle diagram of the fitted modes. The colors 
correspond to the following degrees: blue £ = 0, red £ = 1, 
green £ — 2 and purple £ = 3. The surface of the disks is 
proportional to the heights of the modes. 3-ct error bars are 
overplotted. 



It appears that one of the fitted intervals (v = 855.3 /iHz 
for the £ = mode) has a linewidth somewhat larger than 
what could be expected considering the linewidths of neigh- 
boring orders (see Fig. [TT]). Indeed, the £ = 1 mode of this 
interval (v — 828.2 /iHz) is larger than the other modes in 
the neighborhood. 

Mode heights and amplitudes 



As was mentioned before, we fixed the ratios of heights 
beween modes of different degrees £ in our fit. Because 
HD 49385 is the first CoRoT target for which the £ = 
and £ — 2 ridges are clearly distinct in the echelle diagram, 
we had the opportunity to check the chosen height ratios, 
which are those usually adopted in this type of analysis. We 




E i i i i i =1 

700 800 900 1000 1100 1200 1300 
Frequency (/iHz) 

Fig. 11. Linewidths of the fitted modes with 1-a error bars 
(squares and full line). The linewidths for the alternate 
identification are also given (diamonds and dashed line). 
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performed a new fit releasing the height ratios and found 
Ht=i, n /Ht=o, n = 1.54 ±0.16 
fl"/=a,»/fl<=o, n =0.76 ±0.12 
H e=3 , n /H e=0 , n = 0.20 ± 0.05. 

For I = 1 modes, the theoretical value of the ratio lies 
within the 1-er error bars of the fitted value. For £ = 2 and 
I = 3 modes, the theoretical ratios are at 2 and 3 a of the 
fitted values, respectively. At this stage, this difference is 
too small to be interpreted as an inconsistency with the 
theoretical height ratios. It would be interesting however 
to investigate this matter in future analyses of this type 
on different objects. We noticed that with free height ra- 
tios the results for mode frequencies, linewidths and ampli- 
tudes remained within the 1-tr error bars of the parameters 
obtained with fixed height ratios. 

From the heights and linewidths of the modes we de- 
rived the amplitudes (in ppm), given by A = \ttHT. To 
estimate the error bars on this quantity, we took into ac- 
count the anticorrelation which exists between the deter- 
mination of the height and the linewidth. The mode am- 
plitudes are represented for I — modes as a function of 
frequency in Fig. 1121 They range from 3.6 ppm to 5.4 ppm. 
However, these estimates correspond to instrumental val- 
ues. Taking into account the ins trumental response func- 
tions of CoRoT (as described in iMichel et al. 20091) . they 
can be converted into intrinsic bolometric amplitudes per 
radial mode by multiplying them by a factor 1.032. The ob- 
tained amplitudes are consistent with the prediction of the 
maximum amplitudes using scaling laws (Ai = q = 6.1 ± 0.5, 
see Appendix IA1 . The star HD 49385 is then quite specific 
among the other solar-like pulsators observed with CoRoT 
for which amplitudes were f ound significantly lower 
than the predicted ones (see lAppourchaux et al. 2 008. 
iBarban et al. 20091 iGarcia et al. 20091 IMichel et al. 2 0081 . 

The signal-to-noise (SNR) ratio of each mode computed 
from the fitted height and the background estimated in 
Sect. 15.3. H is given in Table [3] The SNR is found to be very 
small for the fitted i = 3 modes (about 0.25 at the max- 
imum of the signal). However, when releasing the height 
ratios, we find H l= ^ n /H l=0>n = 0.20 ± 0.05. The SNR of 
the I = 3 modes should be closer to 1 at the maximum of 
the signal. 

Splitting and inclination angle estimates 

Most of the fits we performed converge toward a null 
value of the splitting and consequently an undefined 
value of the inclination angle. This can correspond to two 
different cases: 

— The star could be seen almost pole-on, i.e. with an in- 
clination angle close to 0°. In that case the m ^ com- 
ponents of the non-radial modes are not visible, and we 
therefore cannot estimate the splitting. However, the 
measured vsmi of the object and the estimate of the 
stellar radius R allow us to give constraints on the split- 
ting i> s through the relation 

vsvui = 2nRv s sini. (5) 

For a low inclination angle we expect higher values of 
the rotational splitting (e.g. v a > 1.5 /xHz for i lower 
than 10°). 
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Fig. 12. Amplitudes of the fitted modes with l-er error bars 
(squares and full line). The amplitudes for the alternate 
identification also appear (diamonds and dashed line). 

Table 4. Reference results for the parameters of the non- 
radial p modes (symbols are the same as in Table [3]) . 



£ 


7J i) ( 1 1 

^71, £ \t J,11Zj J 


flnr ill „\ 


Hi IB(u A 


1 


777.91 ± 0.24 


* 


3.26 


1 


09c 91 _i_ n /i9 

0Z0.ZI ztz U.^Z 




z.oo 


i 


881.29 ±0.29 




4.44 


i 


935.90 ± 0.23 




6.00 


i 


991.09 ±0.22 




7.03 


i 


1047.79 ± 0.24 




7.62 


i 


1104.68 ±0.22 




5.76 


i 


1161.27 ±0.33 




3.89 


i 


1216.95 ±0.53 




2.01 


2 


794.55 ± 0.52 




1.11 


2 


855.29 ±0.93 




0.97 


2 


905.31 ±0.35 




1.51 


2 


961.47 ±0.49 




2.04 


2 


1017.56 ±0.27 




2.38 


2 


1075.01 ± 0.27 




2.58 


2 


1130.79 ±0.61 




1.95 


2 


1187.55 ±0.32 




1.31 


2 


1246.78 ± 0.84 




0.68 


3 


767.74 ± 0.65 


7.3 cr 


0.11 


3 


821.55 ±1.57 
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928.93 ± 0.57 
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1157.70 ±0.83 
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— We could also explain this result by a low value of the 
splitting compared to the linewidth of the modes, so 
that the fits are unable to separate the components of 
the non-radial modes. Equation[S]shows that a small ro- 
tational splitting corresponds to a high inclination an- 
gle. 

This fitting procedure does not significantly favor any of 
these two hypotheses. However, they lead to almost iden- 
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Table 5. Results of a Lorentzian fit performed on the peaks 
identified as potential mixed modes. 



v (/iHz) 


H (ppmV/iHz) 


r (fiHz) 


a (ppm) 


748.60 ± 0.23 
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tical p mode parameters (all the parameters agree within 
1-a error bars). 

In both cases, the effect of the rotational splitting on 
the observed mode profiles is weak. Indeed, either the an- 
gle is small, and the m ^ components of the multiplets 
have very small amplitudes compared to the m = compo- 
nent, or the splitting is small, and the m components are 
mixed together. We note that this favors the detection of 
lower degree modes and therefore justifies why £ — 3 modes 
seem to be detected in HD 49385, while they were not in 
previously analyzed CoRoT targets. 

5.3.3. Possible mixed modes parameters estimates 

The position of HD 49385 in the HR diagram suggests that 
it either reaches the end of the main sequence or lies shortly 
after it. The star HD 49385 is therefore evolved and its 
acoustic spectrum may contain mixed modes, i.e. modes 
which have a g-mode behavior in the deep interior an d a p- 
mode behavior below the surface (see IScuflaire 19741) . The 
mixed modes are known to experience the so-called avoided 
crossing phenomenon, during which their frequency does 
not follow the asymptoti c approximation a nd they can be 
located out o f the ridges (lUnno et al. 19891). It w as already 
noticed (e.g. iDziembowski Sz Pamvatrrvk h 1991) that the 
detection of these modes could yield a very valuable con- 
straint on the size of the convective core. 

As was noticed in Sect. 15.11 two peaks are detected in 
the power spectrum with a confidence level above 95% and 
are not located within any of the observed ridges in the 
echelle diagram (namely tti and 714). They could be the sig- 
nature of mixed modes. We note that if the peak m, found 
at v ~ 748 /iHz is an £ — 1 mixed mode in avoided crossing, 
it could explain why the I = 1 p mode at v ~ 797.5 /iHz 
seems not to follow the curvature of the t — 1 ridge. We 
fitted these peaks as Lorentzian profiles with free frequen- 
cies, heights and linewidths. Because the mixed modes have 
a g-mode character near the center, their inertia is larger, 
and we expect them to have a longer lifetime and therefore 
a smaller linewidth than regular modes. We obtained the 
results given in Table [5j We notice that the fitted frequen- 
cies have small error bars, but the other parameters are 
ill-determined. 

5.3.4. Other mode identification scenarii 

Can ridge D be the m = +1 component of the £ = 1 modes? 

Still considering the identification preferred in this 
paper (i.e. ridge C corresponding to £ = 1 modes), ridge 
D could also possibly be seen as the m = +1 of the 
1 = 1 ridge, as already mentioned in Sect. 15.11 When 
not including the i — 3 modes in the fit, we indeed find 
solutions with a splitting of about 7 /iHz, corresponding to 
the distance between ridges C and D. 



However, we ran models representative of HD 49385 and 
found that with this rotation frequency the m ^ com- 
ponents of non-radial multiplets are no longer symmetric 
with respect to the m = component. To treat this case 
correctly, we therefore need to take this asymmetry into 
account in our fits. 

For this purpose, we studied this asymmetry on evo- 
lutionary models reproducing the T e fj and L that we ob- 
tained for HD 49385 in Sect. I2,ll and with different angular 
velocities. The mode ls were comp uted with the evolution- 
ary code CESAM2k (jMorel 19971 ) and the eigenfrequencies 
were derived from them with the oscillation code FILOU 
(jSuarez fc Goupil 20081) . They showed that the effects of 
asymmetry should be non-negligible even for values of the 
splitting as low as a few /tHz. We describe the asymmetry 
of the ra/0 components of a mode of degree £ and order 
n by the shift 5 n ^^ m defined as 

c ^n,£,+m "i~ Vn,£,—m, 

0nJ,m = V n ,e,0 ^ ■ ■ (,Oj 

These shifts are represented in Fig.[T2]for a model represen- 
tative of HD 49385 for eigenmodes in the frequency range 
of the detected oscillations (13 ^ n ^ 21) and with a rota- 
tional splitting of 3 /iHz. We observe that the shifts 8 n ,i,m 
are expected to vary approximately linearly with respect to 
frequency. We therefore model the shift 5 n ,e,m as 

S n ,e,m = oti,m(n - n ) + $n ,e,m, (7) 

where no corresponds to the order where the signal is max- 
imal (no = 17). For ra/0 modes in non-radial multiplets 
we now fit eigenfrequencies of the type 

V n ,l,m = ^n,i,0 + mv s ~ Sn,e,m- (8) 

This adds two free parameters to the fit per fitted m ^ 
component. 

The obtained results for the eigenfrequencies, heights, 
linewidths and amplitudes of modes of degrees < I ^ 2 all 
lie within the \-a error bars of those given in Table [3J The 
fit converges toward a splitting of v s = 3.8 ± 0.2 /zHz and 
an angle of i — 20.5 ±2.3°. As expected, ridge D is fitted as 
the m = +1 component of the £ = 1 modes. The obtained 
rotational shift is such that a.£ = i, m= i = 0.30 ± 0.19 /iHz 
per order n, and 5 no ,f=i.m=i = 3.0 ± 0.6 /iHz. However, 
using Eq. [S] and the fitted values of the splitting and the 
angle, we obtain vsini = 11.2 ± 2.8 kms^ 1 . This result 
is not compatible with the firm upper limit of 5 kms -1 
established for v sini in the spectroscopic analysis (see Sect. 

H21). 

Can ridge C correspond to £ = Q modes? 

As was said in Sect. 15.11 for the sake of completeness 
and in the light of previous analyses of p-mode oscillations 
in solar-like pulsators, we decided to also investigate 
the alternate identification. In this case, ridges C and D 
correspond to I = and 1 = 2 modes, respectively. Ridge 
A could not possibly be an £ = 3 ridge, because it has 
much too strong amplitudes. Ridges A and B would then 
be the m = ±1 components of rotationally splitted I = 1 
modes. 

We fitted the modes exactly the same way it was 
done before. The results are given in Tables H2 and [71 
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Table 6. Reference results for the parameters of the radial p modes for the alternate identification (symbols are the 
same as in Table [3]). 
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Fig. 13. Frequency shift & n ,i,m (defined in Eq.[5]) of the m ^ 
components of non-radial multiplets for an evolutionary 
model representative of HD 49385 in the frequency range 
of the detected oscillations and with a splitting of 3 /iHz. 
Crosses stand for I = 1 modes and plus signs for I = 2 
modes. A linear regression of the variations of the shifts 
with respect to the frequency is overplotted (dotted lines). 



and an echelle diagram of the fitted modes is shown in 
Fig. HH The linewidths, heights and amplitudes are very 
comparable with those obtained in the first identification. 
We obtain a splitting of v$ = 2.2 ±0.1 /iHz and an angle 
of i = 70.6 ± 5.8°, which is consistent with the idea that 
ridges A and B are indeed identified as (£ = l,m = ±1) 
modes because they are separated by about twice the fitted 
splitting. Using Eq. [5] as before, these fitted values of the 
splitting and angle give an estimate of vs'mi = 17.7 ± 2.8 
kms -1 , in complete disagreement with the measured 
(wsinz) moasurcd = 2.9t{ r kms -1 . 

The alternate identifications we considered lead to ma- 
jor disagreements with the measured v sini and were there- 
fore rejected. This confirms the identification we adopted 
in Sect. 15.21 to perform the global fit of the p-mode param- 
eters. 



10 20 30 40 50 

Frequency modulo ii/ (/iHz) 

Fig. 14. Echelle diagram for the fitted modes of the alter- 
nate identification. Same symbols as in Fig. [TU] 

Table 7. Reference results for the parameters of the non- 
radial p modes for the alternate identification (symbols 
are the same as in Table [3]) . 
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6. Conclusion 

The star HD 49385 was characterized from both CoRoT 
seismic observations and NARVAL spectroscopic data. The 
atmospheric parameters of the star were derived and the 
137-day-long photometric time series was analyzed. 

A clear series of peaks associated to p modes was de- 
tected in the power spectrum around 1 mHz. Up to now, 
the CoRoT solar-like pulsators presented some ambiguities 
in the identification of ridges corresponding to different de- 
grees I. Here three very clear ridges appear in the echelle 
diagram, which were readily identified as £ — 0, 1 and 2 
modes. The £ — 2 ridge appears clearly distinct from the 
£ = ridge in the echelle diagram. Furthermore, two peaks, 
part of a fainter ridge, were found to be significant and com- 
patible with the characteristics we expect for £ = 3 modes 
in HD 49385. This probably constitutes the first photomet- 
ric detection of I — 3 modes in a solar-like pulsating star 
(other than the Sun). 

We performed a global fit over nine radial orders with 
degrees ranging from I — to £ — 3 modes (simulta- 
neous fit of 36 individual modes). We obtained precise 
estimates of the mode frequencies (uncertainty of about 
0.2/LtHz for radial modes at the maximum of the signal). 
The obtained value of the maximum amplitude for £ = 
modes (5.6 ± 0.8 ppm) is consist ent with the estimate de- 
duced from lSamadi et al (2007\i). unlike the other CoRoT 
solar-like pulsators (e.g. lAppourchaux et al. 2 008). 

We found no evidence of a rotational splitting of the 
modes. This can be explained either by a small inclination 
angle or by a low rotational velocity (inducing a small ro- 
tational splitting compared to the linewidth of the modes) . 
In passing we stress that even for very slow rotators (rota- 
tional splitting as low as a few /iHz) the m ^ components 
of non-radial multiplcts are not expected to be symmetri- 
cal with respect to the m = component. We proposed a 
simple way of treating this assymetry in future analyses of 
solar-like pulsators. 

The p modes of HD 49385 were found to have lifetimes 
ranging from about one day to two days, i.e. somewhat 
shorter than the mode lifetimes in the Sun, but significantly 
larger than those of the previously observed CoRoT solar- 
like targets. This explains why the spectrum of HD 49385 is 
much clearer than for previous CoRoT pulsators, for which 
the large mode linewidths made it harder to separate the 
I = 2 ridge from the £ — ridge. The results obtained for 
HD 49385 confirm that the linewidths of the modes in G- 
type pulsators are smaller than those of F-type pulsators, 
making their analysis easier. 

The very high quality of the spectrum also enabled us to 
detect significant peaks outside the identified ridges. These 
peaks were found to be compatible with mixed modes, 
whose presence can be expected in the spectrum of evolved 
objects such as HD 49385. The existence of mixed modes 
in avoided crossing can explain some specific behaviors we 
observe in the low- frequency eigenmodes. In particular we 
found that the £ = 1 ridge is distorted compared to the 
£ — ridge at low frequency. iDeheuvels fc Michel f2 009) 
showed that this type of distortion could be associated with 
a low-degree mixed mode in avoided crossing. The identi- 
fication of the mode 7Ti as an £ = 1 mixed mode would be 
consistent with the observed pattern. Other features like 
the £ = 2 mode found to overlap the £ = mode (around 
855 /iHz) might also result from avoided crossing phenom- 



ena. This needs to be further investigated in the seismic 
interpretation of this star. 
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Appendix A: Estimate of the expected amplitudes 
for the oscillations in HD 49385 

Scaling laws give an estimate of the amplitude of 
1 = modes at th e maximum of the signal. 
iKieldsen fe Bedding (19 95) established the following rela- 
tion between the luminosity amplitude A of the oscillation 
and the velocity amplitude v osc : 



(A.l) 



We assume that the frequency at the maximum of the signal 
scales as the acoustic cutoff frequency v c , i.e. 
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Inserting Eq. IA.2I and IA.3I into Eq. IA.1I we obtain 
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With the spectroscopic measure of the temperature T c g = 
6095 ± 50 K and our estimate i/ max = 1013 ± 3 /iHz we 
obtain the expected amplitude for £ = modes at the signal 
maximum Ai = q — 6.1 ± 0.5 ppm. 



Among the four studied peaks, two are assumed to be 
£ = 3 modes around the maximum of the signal (it2 and ^3). 
To obtain an expression for P(x\Hi), we need an estimate 
of the expected height and linewidth of £ = 3 modes. 

An estimate of the expected mode linewidths of 
HD 49385 was computed with the non-adiabatic pulsation 
code MAD combined with a ID model reproducing the ef- 
fective temperature, the lumin osity and the log q of the ob - 
ject in the same way as done in lAppourchaux et al. (2 009). 
We obtain around the maximum of the signal a linewidth 
between 1 /iHz and 2.5 /iHz. We adopt a uniform prior 
for the linewidth, taking into account an uncertainty fac- 
tor in the theoretical model: we assume that the maximum 
linewidth is twice as large as that given from the models 
(r e [0,5] /iHz). 

The height of the £ = 3 modes can be deduced from the 
amplitude of £ — modes through the relation 
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for an angle of incl ination i. The visibility factors r^ jTO (i) 
are computed from iGizon fc Solanki ^200^ . The expected 
amplitude of the £ = modes around the maximum of 
the signal is computed in Appendix [X] We impose a uni- 
form prior for the amplitude of the £ = modes (^4 g 
[^4min, Anax], where A m i n and A max are the theoretical val- 
ues of Ai = q at ±3er) and for the angle of inclination (i £ 

[<W2]). 

Based on lAppourchaux et al. (20091 ) we obtain the fol- 
lowing expression for the probability of the observed data 
given the signal 
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where 77 = r max (A max — A m \ n ) tt/2 is the norm alization fac- 
tor. A and v are given in lAppourchaux (20 04*1 and are func- 
tions of the mode height H and linewidth V . 

Inserting Eq. IB. 21 and IB.4I into Eq. IB. II we obtain 
for 7T2 and 773 a posterior probability of 1.4% and 1.1%, 
respectively. 



Appendix B: Computation of the posterior 
probability 

By applying the Baye s theorem and assuming an equiprob- 
ability of Ho and Hi , lAppourchaux et al. (2009T ) obtained 
the following expression for the posterior probability of the 
H hypothesis given the observed data: 



P(H \x) 
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The expression of P(x\Hp) for a sp ectrum binned over n 



ine expression ot r(x\n^) tor a SP ' 
bins is given in lAppourchaux (20041) 

P(x\H ) 
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where 7(71) is the Gamma function. To solve Eq. IB. II we 
need to define the alternate hypothesis Hi and to make 
prior assumptions on the expected signal. 



Two peaks are assumed to be the signature of mixed 
modes (7Ti and n^). For this type of modes, a part of the 
energy is located in the g-mode cavity and the inertia is 
larger. We therefore expect a smaller linewidth for these 
modes. But we can have no clue of their degree £. Because 
both their amplitudes and their linewidths greatly depend 
on £, we use very conservative priors: uniform priors on the 
amplitude of the mode (A G [0, A max ] ppm) and on the 
linewidth (r e[0,5] /iHz). 

We obtain posterior probabilities of 0.001% and 8.7% 
for 7Ti and 7t4, respectively. 



